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Abstract 


The paper discusses the two-loop (NNLO) electroweak radiative 
corrections to the parity-violating Mpller scattering asymmetry in¬ 
duced by insertions to boxes of electron and neutrino mass operators 
(fermion self-energies), vertex functions and boson self-energies. The 
results will be relevant to the ultra-precise 11 GeV MOLLER exper¬ 
iment planned at the Jefferson Laboratory, which will measure the 
weak charge of the electron and search for new physics. The nu¬ 
merical estimations for the NNLO contribution to the cross section 
asymmetry are presented. 

PACS: 12.15.Lk, 12.20.Ds, 13.40.Em. 


1 Introduction 

The Mpller scattering [1] with polarized electrons has attracted active interest 
from both experimental and theoretical standpoints for several reasons. It has 
allowed the high-precision determination of the electron-beam polarization at 
SLC |2], SLAG [2111], JLab |5] and MIT-Bates [6] (and as a future prospect 
— the ILC m)- The polarized Mpller scattering can be an excellent tool 
in measuring parity-violating weak interaction asymmetries [8]. The hrst 
observation of Parity Violation (PV) in the Mpller scattering was made by 
the E-158 experiment at SLAG laiiniin], which studied scattering of 45- 
to 48-GeV polarized electrons on the unpolarized electrons of a hydrogen 
target. It results at = —t = 0.026 GeV^ for the observable parity-violating 
asymmetry Apv = (1.31±0.14 (stat.)±0.10 (syst.)) x 10“'^ [12] which allowed 
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one of the most important parameters in the Standard Model (SM) ~ the sine 
of the Weinberg angle sin 6 ^ 14 ^ - to be determined with accuracy of 0.5 % 

The MOLLER (Measurement Of a Lepton Lepton Electroweak Reaction) 
experiment planned at the Jefferson Lab aims to measure the parity-violating 
asymmetry in the scattering of 11 GeV longitudinally-polarized electrons 
from the atomic electrons in a liquid hydrogen target with a combined stati¬ 
stical and systematic uncertainty of 2 % [a da dsi [IT]. At such precision, 
any inconsistency with the Standard Model (SM) predictions will clearly 
signal the new physics. However, a comprehensive analysis of radiative cor¬ 
rections is needed before any conclusions can be made. Since MOLLER’s 
stated precision goal is signihcantly more ambitious than that of its prede¬ 
cessor E-158, theoretical input for this measurement must include not only 
a full treatment of one-loop (next-to-leading order, NLO) electroweak ra¬ 
diative corrections but also two-loop corrections (next-to-next-leading order, 
NNLO). 

The signihcant theoretical effort has been dedicated to one-loop radia¬ 
tive corrections already. A short review of the references on that topic is 
done in [HI |19| , where we calculated a full set of the one-loop electroweak 
corrections (EWC) both numerically with no simplihcations using computer 
algebra packages and by-hand in a compact form analytically free from non¬ 
physical parameters, and found the total relative correction to the observ¬ 
able asymmetry to be close to —70%. It is possible that a large theoretical 
uncertainty in the prediction for the asymmetry may come from two-loop 
corrections. One way to hnd some indication of the size of higher-order con¬ 
tributions is to compare results that are expressed in terms of quantities 
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related to different renormalization schemes. In (SU], we provided a tuned 
comparison between the result obtained with different renormalization con¬ 
ditions, hrst within one scheme then between two schemes. Our calculations 
in the on-shell and Constrained Differential Renormalization schemes show 
the difference of about 11%, which is comparable with the difference of 10% 
between MS [21] and the on-shell scheme [22]. It is also worth noting that 
although two-loop corrections to the cross section may seem to be small, it is 
much harder to estimate their scale and behavior for such a complicated ob¬ 
servable as the parity-violating asymmetry to be measured by the MOLLER 
experiment. 

The two-loop EWC to the Born cross section (~ can be di¬ 

vided onto two classes: Q-part induced by quadratic one-loop amplitudes 
~ M-iM-i, and T-part - the interference of Born and two-loop amplitudes 
~ 2 Re (here index i in the amplitude Ait corresponds to the order 

of perturbation theory). The Q-part was calculated exactly in [23] (using 
Feynman-t’Hooft gauge and the on-shell renormalization), where we show 
that the Q-part is much higher than the planned experimental uncertainty 
of MOLLER, i.e. the two-loop EWC are larger than was assumed in the past. 
The large size of the Q-part demands detailed and consistent treatment of 
T-part, but this formidable task will require several stages. Our hrst step 
was to calculate the gauge-invariant double boxes [21]. In this paper we do 
the next step - we consider the EWC arising from the contribution of a wide 
class of the gauge-invariant Feynman amplitudes of the box type with one- 
loop insertions: fermion mass operators [or Fermion Self-Energies in Boxes 
(FSEB)], vertex functions [or Vertices in Boxes (VB)], and polarization of 
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vacuum for bosons [or Boson Self-Energies in Boxes (BSEB)]. 

The paper is organized as follows. We define the basic notations in Sect. [2] 
and present FSEB, VB, and BSEB in Sect. [3l In Sect. 01 we provide the nu¬ 
merical results for asymmetry for the kinematics conditions of the MOLLER 
experiment and discuss work still to be done in the future. In Appendix 
the mass operators of electron and neutrino are presented. In Appendix [HI 
we show the result for one-loop corrections to vertex functions for the case 
when only one fermion is on the mass shell. In Appendix [Ul we consider the 
polarization of vacuum for the virtual photon, Z- and lE-boson. The details 
of calculation of ultraviolet cut-off loop momenta integrals can be found in 
Appendix [Dl 

2 Basic notations 

We consider the process of electron-electron elastic scattering, i.e. Mpller 
process: 

e-(pi, Ai) -|- e_(p2, A2) ^ e_(p3, A3) -|- e_(p4, A4), ( 1 ) 

where Aj {i = 1,4) are the chiral states of initial and hnal electrons. The 
kinematical invariants were defined in the standard way: 

s = {pi+ P2)^ t = {pi- P3)^ M = (pi - P4f. (2) 

In the MOLLER experiment, the expected beam energy is Ebeam = 11 GeV, 
that is s = 2mEbeam ~ 0.01124 GeV^, where m is the electron mass {pj = 
m^). For the central region of MOLLER (at 6 ~ 90° in center-of-mass system 
of initial electrons), —t ^ —u ~ s/2 thus we can use an approximation 
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that s, \t\, |m| 3> Also, as for MOLLER kinematics in central region 
s, |t|, |«| -C we neglect in following the terms of order 0{s/mz,w)- 
We consider the process ([T]) in terms of chiral amplitudes Ai^, where 
A = {A 1 A 2 A 3 A 4 } is the chiral state of initial and hnal electrons. The PV 
asymmetry to be measured by MOLLER is then dehned as 

,s^ + t^ + u^ 


A = 


\M -1" - \M- 


= 2 ( 8 ™)^ 


In the Born approximation, this asymmetry has a form 




(3) 


sHu 


2m^ Sm 


(4) 


w 


proportional to 


a = 1 — 4s 


w 


(5) 


Let us now recall that sw (cw) is the sine (cosine) of the Weinberg angle 
expressed in terms of the Z- and lE-boson masses according to the Standard 
Model rules: 

Sw = \/l - C^, Civ = miv/mz. (6) 

Thus, the factor a is just a ~ 0.109 and the asymmetry is therefore suppressed 
by both sjvn^ and a. Even at at 6 ^ = 90°, where the Born asymmetry is 
maximal, it is extremely small: 

r(o) = —~ 9.4968 • 10"®. (7) 

We denote the specihc contribution to the asymmetry by the index C, which 
thus can be BSEB, FSEB, VB or IB=BSEB+FSEB+VB for the whole set 
of diagrams Fig. [H respectively. 
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The contribution to the asymmetry (AA)^ and the relative correction 
are dehned as: 


(AA)c = 


nC” _ 
— 


|A^-— 


^ | A ^^|2 


(AA)c \M 


0 I 
c I 


\M 


c 


A(o) \M^ -P-I-Mc 


( 8 ) 

(9) 


The relative correction to observable asymmetry from the contribution 
of type C looks as (see derivation in more details in [T^): 


- 


_ ^(0) 
A(o) 


1 + 5^ ’ 


( 10 ) 


where the relative correction to unpolarized cross section ctqq (we used short 
notation for differential cross section a = da/d{cos6)) is: 

a° 

^00 

For the two-loop effects where 6^ is small, we can use an approximate equa¬ 
tion for relative correction to asymmetry 6^ ~ D^. 



3 Insertion of mass operator, vertex and vac¬ 
uum polarization functions to the box type 
amplitude 

The numerical value of loop momentum squared |/c^| in the box-type ampli¬ 
tudes with the heavy boson exchange is large compared with the square of 
electron mass \k‘^\ ^ since if |/c^| is far from the contribution is 

suppressed with the mass of heavy boson squared in denominator. So we can 
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(a) 


(b) 


Figure 1: (a) Fermion self energies in boxes (FSEB), (b) boxes with vertices 
(VB), (c) boson self energies in boxes (BSEB). In this diagrams all wavy 
lines are assumed to be photons or Z-bosons. We also considered crossed 
box diagrams which can also contain IF-bosons legs. 

use the asymptotic expressions for the one-loop vertex functions as well as 
the mass and vacuum polarization operators. Using the well-known approach 
pS] . ^2E\ which we successfully employed for the box-type chiral amplitudes 
in [23] (see also [2Z]), we can write for the direct ZZ-hox chiral amplitude of 
“++++” type: 


uslf^ia + 75)^1 ^47^(0 -F 75 )(-^) 7 ''(« + 75)^2 





Easily we can get a similar expression for the crossed box and for amplitude 

of “-” type. The quantities g and / in (fT 2 |) coincide with a and b from 

[23], respectively, and are defined as: 


g = UiUJ_P2UJ+Ua, f = U2UJ_piUJ+U3, 


where a;± = (1 ± 75 ) /2 are the chirality projection operators. Let us calcu¬ 
late f-channel amplitude; the u-channel amplitude can be obtained by the 
replacement f 77 u. This interchange will be denoted below as an operator 














The box-type amplitude with the double Z-boson exchange with all the 
possible insertions (i.e. VB, FSEB and BSEB) has a form: 


M 


zz,± 


OO 

. 0:^(1 6s^t C dr . 

* gfml J 


( 13 ) 


0 

where “±” sign corresponds to the chiral amplitudes . The expression 

for the box amplitude with Zy-exchange is similar: 


OO 

A^^7,± _ •2a^(l ± f dr 

^ iAcwSwV gfml J r(l + r) ^ 

0 

At last, for 77 -exchange amplitude we have: 


= 


6 sH 


la 


gfm 


T ( \ 
_2 -^77 (''’)• 


r- 


(14) 


(15) 


In all the above cases, the integration variable is related to the loop momen¬ 
tum as r = —The lower limit of integration z = —t/m| for Ad'*"*' 
is introduced to avoid the double counting for the region of small loop mo¬ 
menta squares —k‘^ < s, where we use the Yennie-Frautchi-Suura approach 
[25] . Finally, the contribution to Ai arises from the box-type Feynman 
diagram with two hF-boson exchange: 


■^ww = 


. o? sH f drw J , N 

—T \ TTh- ^^wwi^w) 

2s^gfmlrJ [l + Tw)^ 


Tw = -k^/ml,. ( 16 ) 


The structure of the quantities Aj in (IT^ . (ITTll . (|T5ll and (1T6|1 corresponds 
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to three types of radiative corrections, FSEB, VB and BSEB, respectively: 
^zz = 2 A^e + 414eZ + 2112 ^, 

-fz 7 = 2A^e + 214eZ + 214e7 + ^ZZ + 14^7, (17) 

*^77 ~ 27Vle d" 4V^e^ 20^^, 

7wrv = 2A4,y + 414j^vk + 2nvKrv- 

Here, we nse the dimensionless qnantities for the prodnct of fermion Green 
fnnction and the trnncated mass operators of electron Me and nentrino My 
(see Appendix EJ: 

1 n 

Me,. = ^Me,,. (18) 

The vertex fnnction with one electron on the mass shell and another 

electron off the mass shell is normalized as 


= -ter Ke7(0) = 0. (19) 


The vertex fnnction is normalized at the point /c^ = m|: 


KeAk") = 


le 


4 :C]yS\Y 

and similarly for ez/lE-vertex fnnction we have: 


r Kez(fc'), Kez(m|) = 0, 


le 


Kw(n = ^=— V,M-(ml-) = o. 


( 20 ) 


( 21 ) 


\/ 2 sw 

The explicit expressions for the vertices Vee'yi Kez and Veuw are given in 
Appendix [Bl 

The dimensionless prodncts of boson Green fnnction with the relevant 
regnlarized polarization operator + B{q^)q^qy are dehned 
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as: 


n, = (,^). 

Hz = 


n“ (0) = |in“ (0) = 0; 

nzz(^z) “ ^j^^zz(^z) “ Qj 


nzz = ■^ng,(4r"), n‘z',(0) = o; 

% d 

^ „ 2_™2 ^wwi^'^)y ^wwi^w) = ’^^wwi'^w) = 0 - ( 22 ) 

The structure B{q‘^)q^q^ does not contribute due the gauge invariance. The 
explicit expression for the “truncated” quantities are given in Appendix O 


4 Numerical results and conclusion 

For the numerical calculations, we use the central kinematical point of the 
MOLLER experiment and a, mw and in accordance with the Particle 
Data Group [29]. The effective quark masses used for the vector boson self¬ 
energy loop contributions are extracted from the shifts in the fine structure 
constant due to hadronic vacuum polarization (^|) = 0.02757 [30] . 

For the mass of Higgs boson, we take mn = 125 GeV. 

The contribution relevant to the observed asymmetry is the interference 
of the two-loop box-type amplitudes with the Born amplitudes The 

contribution to the matrix element squared (i.e. cross section) has the form: 

= 2 (1 + Ptu) . (23) 

In the right-hand side of this equation, we assume that the amplitudes 
are taken in the same chiral state corresponding to the state of left-hand 


11 




side. Note that the intermediate states with bosons and Faddeev-Popov 
ghosts contribute to the mass and vertex operators in the AA chiral 
amplitude. Since the parameter a is very small, we can present the hnal 
result as: 


M 


IB 




" = -H{a) + {H{-a) + Y) = 


da 


a^O 


+ Y, 

(24) 


and thus the relative correction has the form: 


2,,2 


^IB _ __ 

^ 128 (ttq;)^ (s^ + 

We dehne H and Y as: 


— 2 a 


dH{a) 


da 


a^O 


+ F 


WO) 


(25) 


H — Hzz + Hz^ + + Hww + Ffmix, (26) 

Y = Yzz + + Y'i'i + Yww + Wiix, 


where the hrst four terms in both H and Y correspond to the box-type 
amplitudes with ZZ, Z'y, 77 and WW bosons exchanged between electrons, 
and the last term corresponds to the cases with Z or 7 and the mixed boson 
Green function with polarization operator Tiz'y- 

Using the following relations (see, for example, [21] and 1211) 


gf \9f cd) st^u 


1 r t uV _ 2 
gf \9f cd) sH 


(27) 
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we obtain the following numerical results: 


3a^TTil + a)^ n I P ^ 
nzz - - ^7 - 77 -+ ^tu) 


dr 


{l + ry 


Yzz — 


8{cwSwy rn\tu 

0 

+2Uzz + + (1 + afV^)] = 

= -1.653 ■ 10"^=^ (1 + a)^ (-81.36 - 1.1293 (1 + af) 

OO 

Sa^TT . f dr 


[ 2 (AC + (1 + )+ 


/^z^ = 


8{cwSw)‘^ m\tu 
= 3.139- 10"^^; 
12Q;^7r(l + aY 


(1 + Ptu) 


{l+Tf 


[ 2 Mf + A{V^ + = 


{cwSwY rrTztu 


(1 + Ptu) 


dr 


r(l + r) 


[2(M; + (l + a)X^) + 


+ + n,, + 2(1-^ + K" + (1 + af{V^ + VYJ)] = 


= -9.155 ■ 10"^^ (1 + a)^ (-4.30744 - 0.04567 (1 + a)^) ; 


Vzy = 


12Q;^7r 

(civsiv)^ m%tu 


(1 + Ptu) 


dr 


r(l + r) 


[ 2 Mf + 2{V'^ + 1 / 2 " + Vjl)] 


= 6.974- 10"^\ 
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H.. = + (1 + afM-H 


(CwSw) 


' m\t J r2 


+ 2U,, + 4{VZ, + {l + a)X%)] = 

= -3.094 ■ 10"^2(1 + a)2(-2.52038 - 5.04456 ■ 10"®(1 + af), 

OO 

K, = -i^(l + F,„)5r f + 4VZ] = -4.4261 ■ IQ-'', 




m|t J 


Hww — 0 ; 

Sa^TT 

^ww = 


{swY m\tu 
= -3.36 ■ 10"^°. 


(1 + Ptu) 


dr\y 


(1 -h TwY 


[2Mi, + 21lww + P^eWu] = 


(28) 


The “mixed”-type amplitude in two-loop approximation has two different 
contributions (if, T)mix = (if, + (-^i ^)mix- The hrst contribution is 

associated with the two-loop box-type amplitude: 




(cw-sw)^ m^tu 


0 


(1 PTw)' 


1 


X ( {AcwSwY' -^ (1 + - 

V T-w l + T-w 

= -1.10029 ■ 10"®(1 + a)(0.007746 - 0.000340(1 + a) 


= 0 

mix 


R'lzi'T'w) — 




g —3-h 2(3 — 2c^)— ]Lw{tw), 
bTTSvn \ 'Tw, 


Lw{tw) = J dx\og{l + x{l — x)tw)- 
0 


The second contribution arises from the interference of the Born-type am- 
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plitude with the mixed Green function and the box type one-loop amplitude 
with the 77 -exchange: 



3.982-10-‘®(1 +a), 


(29) 


The contributions to the asymmetry from the transition polarization op¬ 
erator Tlz'-i with leptons in the fermion loop are proportional to higher powers 
of a, which is small. The same reasoning is valid for the quark-antiquark state 
contribution. Specihcally, it enters with the factor 


(2/3)(l - (8/3)4-) - (1/3)(1 - (4/3)s?,,-) = a/3. (30) 


The contributions from (G^G^) intermediate states 

are considered in Appendix O 

Finally, we are ready to present hnal numerical value for the relative cor¬ 
rections considered in this paper to the observable cross section asymmetry. 
The one-loop (NLO) corrections [iHl HE] give the biggest contribution, 


^NLO ^ _o.6953. 


(31) 


Several categories of the NNLO contributions (Q-part and double boxes) are 
calculated in [23] and [21] and give the following values: 


^NLO+Q _ _o.6535, ^ T^doubie box _ _0.0101. (32) 


Summing up all the contributions in fl25|) . the numerical result of the class 
of the gauge-invariant Feynman amplitudes considered in this paper (boxes 
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with one-loop insertions of fermion mass operators, vertex functions and 
polarization of vacuum for bosons) is: 

^ = -0.0039. (33) 

As one can see, the relative correction we obtained is much less than the 
expected MOLLER experimental error, but it still a non-negligible contribu¬ 
tion to the MOLLER error budget. Most likely, the entire set of two-loop 
corrections will be smaller than the experimental statistical error, but, in 
the light of the MOLLER success depending so crucially on its precision, the 
two-loop corrections still need to be controlled. 

As the low-energy precision experiment, MOLLER is complementary to 
the LHC efforts and may discover new physics signal that could escape LHC 
detection. However, for the MOLLER experiment to produce meaningful 
physics, the uncertainties in the NNLO EWC must be much smaller than 
the MOLLER statistical error. Clearly, there is a need for the a complete 
study of the two-loop electroweak radiative corrections in order to meet the 
MOLLER precision goals. 
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A Mass operators 

Here, we will define the explicit form of the quantities Mg, My which enter 
to lij from flTTl) . The quantity Mg has the following form: 

Mg = M2 + (a ± 75)'Mf + ujImY • (34) 


The explicit expression for the truncated mass operator in QED was found 
by R. Karpins and N. Kroll in 1950 [3T1 [32] : 


M2 


\2 

[p — m) 


271171 
p + m 


[2(1-P) 


1 2 3p I 

1 -y—-logp) - 


mp V2(l “ P) 


2 -p+di2^1ogd+l + 21ogA 
1 — p J m 


(35) 



It is useful to note that the expression in the square brackets is hnite at 
p ^ 1. In the limit of large ti = with logarithmical accuracy we 

have 


CK 

M2 = M2 (ri) ■ p ^ log (ti) • p, Ti > 1. (36) 

dvr 

This mass operator contribution to the integral in 0151) with logarithmical 
accuracy gives: 


OO 



-t/m| 
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The mass operators induced by additional Z and W bosons have the following 
form: 

1 

0 

1 

Ml, = MY = [(1 — a;) log(l + Tx)dx. (38) 

J 
0 

B Vertices 

The general form of the vertex function is V^{k) = A'y^ + Bk^] the term Bk^ 
inserted in the box-type amplitude gives no contribution due to the gauge 
invariance. The vertex function with one electron on the mass shell and other 
electron off the mass shell — k,k) = —iej^Vee-yi^k'^), normalized as 

14 e 7 ( 0 ) = 0 , has three contributions: 

= Ky + (o ± 75 )VJ, + (39) 

First, let us consider the QED-type contribution with the virtual photon 
intermediate state The standard procedure of joining denominators 

and performing the loop momenta integration leads to 

,24 « f f f M k%h\ III. 

) = — dx ydy nog— + 1 , b = xy, b = 1 - b, 

0 0 

(40) 

D = {m^ — k‘^x{l — x))y‘^ + (1 — — 1/(1 — y){k'^ — ‘^Pik), 

where A is cut-off regularization parameter. Since the sub-set of the diagrams 
considered here is gauge invariant on its own, it was not essential for us to 
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use the dimensional regularization scheme providing gauge invariance, so 
we simply applied the cut-off technique. There is no signihcant numerical 
difference between two schemes in this situation. 

The renormalization procedure consists in subtraction at /c = 0 and leads 
to: 

1 

KI 7 i^e) = j log(l + x{l - x)Te)dx, Te = (41) 

0 

The contribution of this vertex function to the integral in flTSl) has the form: 


dr 


a 


nir 


^ ( 1 - log — ) . 


dvr 


-tlm% 


—t 




(42) 


The other contributions are: 

1 1 


^ - 

27i{4:CwSwy 


dx / ydy log 


1 - 2 / 


bbr 


1 — y + bbr 2(1 — y + bbr) J ’ 


yW ^ ^ 

dvrs?, 


w 


0 0 


(43) 


0 0 

Vertex function = —iG 7 ^(a± 75 )\ 4 ez, G = e/{AswCw) has four different 


contributions: 


Veez = + (a ± 75 ) + a;_ V" + (44) 
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and is normalized as Veezik"^ = ^“z) = 0. These contributions are 


a 


r" = -j^iogT, 


= 


a 


{AcwS\y)'^ 2tt 


X log 


dx J ydyx 
0 

bbr 


bb 


I — y — bb 

1 — y + bbr 2(1 — y + bbr) 2(1 — y — bb) 
1 1 


= 


ac 


w 


dTT 


X ( 3 log 


J dx J ydyx 
0 0 

yc^ + Tbb Tb{b — b) 


b{b - b) 

ycl^ - bb 2{yc^ + rbb) 2{yc^ - bb) 
1 1 


1 / 2 *^ = 


a 


27rs^ 


dx / ydyx 


X log 


(1 - y)cl^ - bb 


rbb 


bb 


( 45 ) 


(1-9)4'+ ’■55 2(9C^,+ tW)) 2(ycly-hb)^ 

And finally, the vertex function as well contains three 

contributions: 


14 ^^ = (46) 


and is normalized as Veuwi'!' = ~c^) = 0 and So the contribu- 
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tions are: 


yZW ^ 


a 


4vrs^ 


dx / ydyx 


X 


„, ycLx + Tbb rbib — b) 

-3 log ^ ^ ^ 


yjW ^ 

47r 


yax - 
1 1 


cl^bb 


c^b{b - b) 
2{yax; + Tbb) ' 2{yax, - c^bb) 


+ 


dx dy X 


.0 0 

31 + r{b-b) 

bcw 


m 


l + ^log^ + 




1 m? 


2(c^ + r6)y ' 4 
where ax = x + (1 — x)c^. Note that the term containing log(m^/A^) in 
expression can be omitted as it will be absorbed by the similar terms in 
two-loop contributions after applying the Yennie-Frautschi-Suura regulari¬ 
zation (see [33] for details). 


C Polarization operators 

While considering the vacuum polarization operators of photon, Z- and W- 
boson at one loop, one should recall that the regularization implies the double 
subtraction procedure. The “truncated” operators imply including only the 
vertices of interaction of bosons with the fermion loop. From now on, we will 
omit index “tr”. The general form of the polarization operator is: 

= gpiv^{.q^) + q^,quB{q^). (47) 

We only need to consider a part of polarization tensor proportional to g^y. 
The reason is the gauge invariance of the whole set of the double-box amplitu¬ 
des, which leads to a zero contribution for terms proportional to q^^q^ tensor. 
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Let’s define IT'’' as: 


( 48 ) 

It has five types of contributions, corresponding to the intermediate state 
of lepton-antilepton pairs, quark-antiquark pairs, W^W~ and the charged 
ghost state G^G^\ 

m = + + + + ( 49 ) 


The contribution of leptons and quarks are associated with the quadratic 
divergent integral over the loop momentum: 


dk 


-Sp (fc + - g + m)7i, 


(50) 


A J (k^ — rn?) {{k — g)^ — rn?) 

Using the set of divergent integrals (see Appendix ED and performing the 
regularization procedure, we include the contribution of leptons and quarks 
as 

a 


n‘ + m = 


TTT 


(5^G(T,ff') + 3 ^ Q2G(r,a’)j, 

\l=e,fi,T q=u,d,s,--- / 


(51) 


where 


G{r, (t) = hT- 2a)L{-) + ^r, L{z) = [ dx log(l + x(l - x)z), 
3 cr 9 ./ 


/ = !!u 

C 2 ) 


T = 


q 


m 


2 ’ 


The factor 3 takes into account the number of quark colours. The last three 
contributions in (1491) are 


Y^ww ^ 




a 


127 rr V 6 


^r+(5r-c^)L (^ ) ) , (52) 
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the known result for Feynman-t’Hooft gauge used in [SUES]. 

The polarization operator for Z-boson has seven types of contributions: 

+ n| + , (53) 

were we used the dehnition 

(54) 

q UL^ 

The contribution of lepton 11^^, quark n| and the neutrino fl^ loops can be 
calculated in the non-renormalized approach: 

^ (^q'^ log ^ + 0(g^)^ gi,u- (55) 

The renormalization of R{t) for any contribution to the polarization operator 
of Z-boson consist of the following replacement: 

R{t) R{t) - R{-1) - (r + l)i?'(-l). (56) 


In particular, for example: 


Q 

g^log-^ ^ m|F(r), F(r) = r logr - 1 - r. 




(67) 


Keeping in mind that there are three generations of charged leptons, neutri¬ 
nos, and quarks, we obtain: 


n 


l+q+v _ O: 

12tt 1-\- t 


3 + 


+ 


4:{swCw)^ 2(siyCiy)^ 


1 _ 2^2 + 

-L ^ g 


(58) 
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The contribution of W~^W pair in the intermediate state to the Z-boson 
polarization operator looks like: 


n 


w+w- 


ac 


w 


Stts^ 1 + r 


X 


X 


/19 16\ / r \ /19 16\ 17, 

( —^ ^ ^ ) Cl + —(r + 1)C2 


V 6 


'W 


Ci = L{ —^ 




-0.248, 


6 


C 2 = 


The contribution of the charged ghosts is: 


2 

x{l — x) 


1 — x{l — x)/c 


dx ^ 0.226. 


yt^w^w _ 

‘-‘-z — 


X 


- 

‘-‘-Z — 


a{l-2sl,f 1 
4:7l{cwSwy 1 + T 


12 ^ + s'"' 
1 


X 


T 




“ + 1 ) 




- gCiy ) C2 


27r 1 + r 


L ( ] — Cl-^C2(l -f t) 

CwJ c^ 


And, finally, the contribution from the state with ghosts Gi ,2 is: 
a 1 


n 


G 1 G 2 


-[r(A(r)-A(-l)) + (r + l)A'(-l)], (59) 


4Tl(cwSwy 1 + T 
with explicit form of A{t) given in Appendix [Dl 

The polarization operator for IT-boson has contributions from the loop 
Feynman diagrams with (z/, e), (d + s){u + c), (IF, Z), (IF, 7 ) and the states 
with ghosts. Defining the dimensionless combination: 

-.2 


n'^ = 


nf,(g2) = n'^(riy), Tiy = 


tW 




Q 


m 


2 ’ 
w 


( 60 ) 


we write 




Gw,l 


( 61 ) 
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From now on, when considering the dehnite contributions to we imply 
that r —)■ Tw- Let us hrst consider the contributions from fermions. For the 
state with a charged lepton and the corresponding antineutrino we obtain: 


En 


Ivi _ 


a 


“ ^4x4,1 + 


( 62 ) 


with function F given in flSTl) . Factor 3 corresponds to the number of lepton 
generations. The contribution of quark states is: 

1 


— 4- 

q=u,d,s,c 


a 


247rs^ 1 + r 




(63) 


where factor 4 corresponds to the number of pairs (d + s)(m + c). The for the 
WZ state we have: 


T\WZ _ 
— 


ac 


w 


4s^(l + tiy ) 


[vI/(riy)-vl/(-l)-(l + r^)vI/'(-l)], 


(64) 


4 /( 4 ) = (az — 1 - y j [ log I X + ^ 2 ^ + “ x)z\ dx- 

V J J \ J 


- ( ^^ + ^ ) / log(l + x(l - x)z)dx+ 

1 1 

+ j j y (^ 1 / + (1 - y) + y^^ ~ y)^ ^y^ 


0 0 

T(-l) = 0.226, T'(-l) =-1.26. 


Now we consider the intermediate states {W^Gz) and {Gw^Gz)- For the 
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insertion to the box amplitude we have: 


T-fWGz I tjGw^Gz 


a 1 
27r 1 + r 


dx log 


a; + (1 — a;)(l + tx)cI 
x + {l- xyci, 




with A{t) taken with 7 = 1/c^. 
For the last two terms we have: 


tM^.7 I tt'^W.7 _ 


tt|/V,7 I ttLivv 


47r(r + 1 ) 


-4<5(r) + ) > 


Q{t) = t [ dx log(l + rx) + 3 + (1 + r) f 1 + ^ log ^ 


1 ^ 

/?(r) = - 4 - —+ 11 + 6 


r 


1 + r 


log(l + r) - 20 - 27(l + r), 


Note that the term log {m?/\^) in the expression for Q{t) is compensated by 
the corresponding contributions from the two-box amplitudes. 

Let us now consider the contributions to the transition polarization = 
and dehne the dimensionless function 




As shown above, the fermions contribution is proportional to and can be 
omitted. The contributions of (W^W~), {W^G^), {G^G^) to 11^'^ are, 


respectively: 


. acw f 19 16 

^Sttsw V 6 3tu 


L{tw), -i 


. acw f 1 


Sttsw V 6 3tu 


L{tw)-, 


27rsw'+v 


L{tw)- 
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Thus, the total is: 


f-3 + —(6 - 44 ,)) L(r„). (67) 

StTSw \ ) 


D Loop integrals and regularization 


To calculate loop integrals, we perform the Wick rotation of the loop mo¬ 
mentum k (fco —^ *^ 4 , k'^ = < 0). In order to regularize ultra-violet 

divergence, we introduce the cut-off parameter A so < A^, and all of the 
kinematical invariants much less (i.e. A^ S> \PiPj\)- The final result will be 
independent of A after the renormalization procedure. Let us now list all the 
integrals we need: 


k‘^dk 


{k^-Df 

dk 

{k^-Df 
{efdk 
(F - DY 


log^ 

1 

3 

f dk 

= log^ 

1 

2’ 

J [k^-Df 

r dk 

2D' 

\og- 

11 

J {k^-ny 

f k‘^dk 

6 T)2’ 

1 

6 ’ 

to 

1 

3D' 


( 68 ) 


Here, we use the notation dk = d^k/ (fvr^) = k'^dk'^, where kE is the Euclidean 
4-vector (i.e. k^, = k\+k 2 +k‘l+k‘l > 0) and omit the terms of order 0{D/ 

We also use the consequence of the integrand symmetry: 


/ 


/ k^dk = 0 


(69) 


for any function /(A:^). The standard procedure of shifting variable in loop 
integrals [52] leads to: 


dk A^ 

{{k-hf-df ^ T “ ’ 


k^dk 


{{k-hf-dy 
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Let us consider the divergent integrals with A = k‘^ — m? and B = {q — k)"^ — 
2 

m . 


dk 

~AB 


— La — 1 — L, 


k^dk 1 [ 3 

= - - - L ) , 


AB 2 


k^kydk 


AB 


A2 


q m 


72 


- 1 - + ;;::^- :r + \ m - — ] L/^ + - { - - m ] L > + 


6 


1 fq 


3 V 4 


where 


+ qi^qu { ^La - ^ ^ (m^ - g2) L }> , 


f q^ 

LA = log—, L = L{t)= I dx\og{l +x{l- x)t), t = - 

j 


(70) 


By contracting indices in the tensor integral flTU]) . we obtain: 

/ k'^dk , 9 9 9 ^ 9 ^ 

—— = —A^-+ 2m^LA — rn?L. 

AB 2 


(71) 


According the renormalization procedure, we can omit terms having the form 
aq^ + bm? and (cg^ + drn?)L\. 

Let us consider now the general integral of the form 


T _ f kf^kydk 

J {k"^ — m\){{k — qY — rri^) 


(72) 


Now, let us use the following algebraic identity: 


2qk — q^ 


+ 


{2qk — q 


2\2 


{q — ky — ml k'^ — ml {k‘^ — mlY {k‘^ — mlY{{q — ky — ml) 

(73) 

Due to our renormalization convention, we can omit the hrst and the second 
terms in the right-hand side of this equation so the integral reads as: 

k^kty{2qk — q^^dk 


d my 




(/c^ — ml){k‘^ — ml)‘^yk — g)^ — ml) 


(74) 
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First, we combine the factors {k‘^ — m‘1) and in the denominator 

using the Feynman trick: 


1 ^ r (1 — x)dx 

aP'b j (a(l — x) + bx)^ 


(75) 


and obtain 


1 


1 7 (1 — x)dx 

j (/c^ —M^)^’ 


Ml = (1 — x)rnl + xrnl. 


0 

Next, we join the resulting expression with the factor {{k — qY — fnl) with 
the similar Feynman identity: 


f Y-yfdy 

c^d J (c(l -y) + dyf 


(76) 


and, hnally, get: 


1 ^ 7 (1 — yYdy 

{k^-MlY{{q-kY-ml)=^J {{k-yqY-midY^ 

where 


(77) 


d = Tiy(l -y) + ^ = x{l - y) + 7 | 9 + (1 -x)(l -»)], 

9 2 

q mi 

Tl = - 2 ’ 7 = —• 

mf mf 

Thus, we have the logarithmically-divergent loop momentum integral, which 
allows the operation of the loop momentum shifting k = k + qy. After that, 
we can use the loop integrals from the beginning of this Appendix. Now, we 
have: 


= A(ri,7) <f9ixu + 0{q^qY, 

1 1 

^(^,7) j j dyY--x){l-yf (^logd- 

0 0 
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therefore the renormalization procedure for this integral has the form: 

riA(ri, 7 ) ^ - yl(-l, 7 )) + (1 + ri)A'(-l, 7 ), (79) 

where A(— 1 , 7 ) ks —0.0896 and A'(— 1 , 7 ) ^ 0.00654 for 7 = m\lm\ = 
1.879. 
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